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We study the effect of commensurability (integer filling 
factor) on the superfluid (SF) - Bose-glass (BG) transition 
in a one- dimensional disordered system in the limit of weak 
disorder, when the effect is most pronounced and, on the other 
hand, may be traced via the renormalization-group analysis. 
The equation for the SF-BG phase boundary demonstrating 
the effect of disorder-stimulated superfluidity implies that the 
strength of disorder sufficient to restore superfluidity from 
Mott insulator (MI) is much larger than that enough to turn 
MI into BG. Thus we provide an explicit proof of the fact that 
at arbitrarily small disorder the SF and MI phases are always 
separated by BG. 

PACS numbers: 05.30. Jp, 64.60Ak, 67.40.-w 

The problem of superfluid - insulator transitions in 
commensurate and disordered bosonic systems has at- 
tracted a lot of interest in the past decade (for an intro- 
duction see Refs. A special part in the analysis of 
the problem is played by one-dimensional (ID) systems, 
where along with the exact finite-size numerical treat- 
ment (quantum Monte Carlo j|,|5| and exact diagonal- 
ization |||7|) there exists macroscopic renormalization- 
group (RG) description [Q. [Recently it was demon- 
strated that the two are in an excellent agreement with 
each other at least in the case of SF-MI transition ||.] 
Note also a remarkable success of density-matrix RG 
study of ID bosonic Hubbard model ||. 

Despite the fact that now the understanding of super- 
fluid - insulator transitions is close to being complete, 
especially in ID, there is one point, however, which re- 
mains somewhat unclear up to now. That is the ques- 
tion of superfluid - insulator transition in a weakly dis- 
ordered commensurate system with integer filling factor. 
From general considerations one would expect here that 
SF and MI phases are always separated by BG phase 
fsLfloJ, while a number of approximate and exact (but 
finite-size!) methods suggest that small disorder is not 
relevant and there occurs a direct transition from SF to 
MI JTT| , p^| j9|] . In this paper we consider this problem in ID 
taking advantage of the asymptotically exact RG treat- 
ment. 

Our main result is the demonstration of the existence 
of the BG phase at arbitrarily weak disorder. We provide 
the proof of this fact in a 'direct' and in an 'enhanced' 
forms. By 'direct' form we mean tracing the renormal- 
ization picture (in the superfluid phase, where the RG 
equations are asymptotically exact) with the observation 
that the commensurate scenario, taking place at smaller 



scales of distance, is destroyed at larger distances by suf- 
ficiently strong disorder, in favor of the universal disor- 
dered scenario. [Incidentally, this accounts for the effect 
of disorder-stimulated superfluidity.] So that commensu- 
rability turns out to be irrelevant to any second-order- 
type transition arising with decreasing of the strength of 
disorder. The 'enhanced' proof does not rely on any as- 
sumption concerning the type of the transition leading 
to the destruction of superfluidity. We simply take the 
phase boundary for superfluid - insulator transition fol- 
lowing from the 'direct' treatment as a lower boundary 
for the strength of disorder sufficient to support superflu- 
idity, without specifying the phase adjacent to SF. Then 
we compare the result with that for the MI phase bound- 
ary, following from the exact relation between the mini- 
mal amount of disorder destroying MI and the insulating 
gap of the pure system The comparison shows that, 
in relative units, the separation between the boundaries 
is ever increasing with vanishing disorder, thus revealing 
the presence of BG between MI and SF. 

The RG descriptions of the superfluid-to-insulator 
transitions in ID were initially constructed in terms of 
the effective action for the long-wave density deviations 

f, employing Haldane's representation of ID bosons 
For our purposes, however, it is more convenient 
to use a dual treatment in terms of the effective action 
for the phase field jlj] . This language directly takes into 
account effects of the collective backscattering off either 
commensurate or disordered potential, the corresponding 
instantons being just vortices in the phase field in (1+1) 
dimensions. As a result, the descriptions for the com- 
mensurate and disordered cases turn out to be very close 
to each other. The only difference is that in the dis- 
ordered case one deals with 'vertical' vortex-antivortex 
pairs (the spatial coordinates of vortex, x\, and antivor- 
tex, X2 , coincide up to some microscopic length) , because 
contributions of non-vertical pairs are averaged out due 
to the random phase associated with the pair: 



l(xi,x 2 ) 



-2tt 



no(x) dx 



(1) 



where no(x) is the expectation of the density at the point 
x p3. In the weak-disorder limit no(x) is related to the 
disordered potential, W{x), perturbatively: 



n (x) =n ( -°\x)+ I dx' S{x',x)W{x') 



(2) 



where n^(x) and S(x',x) are the density expectation 
and static susceptibility for the pure system. 
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S(x',x) = i dt ( [n(x',t),n(x,0)] ) , (3) 



n is the density operator. Substituting (|^) into (Q) and 
taking into account that in the integer-filling case the 
contribution from n^°'(x) leads only to some irrelevant 
renormalizations Ej], we get 



l(xi,x 2 ) = -2tt / dx Q(x)W(x) , (4) 



where 



Q(x) = I dx S{x,x) 



(5) 



is a periodic (commensurate) function. [In Eq.([l|) | 
X2 — x\ I is meant to be much greater than the micro- 
scopic static correlation radius characterizing the decay 
of S(x', x) as a function of | x — x' |.] 

For defmiteness, we take W{x) to be bounded, 

- W/2 < W(x) < W/2 , (6) 

and symmetric with respect to the change of the sign. 
Introducing the radius of a vortex pair 



(7) 



2/i (2/2) being the coordinate of the vortex (antivortex) 
along the imaginary-time axis (y — ct, c is the sound 
velocity, r is the imaginary time), we see from Eq.(Q) 
that the mean-root-square value of the phase associated 
with a pair of the radius R is 



1r 



WVR 



(8) 



[From now on we treat W and R as dimensionless quan- 
tities, assuming them to be measured in units of char- 
acteristic interparticle energy and distance, respectively.] 
While jr <C 1, the disorder is irrelevant and the renor- 
malization picture follows the scenario of the pure com- 
mensurate case. At the scale 



R* 



1/W 2 



(9) 



corresponding to 7^ ~ 1, the renormalization picture 
crosses over to the disordered scenario which becomes 
well-developed at R ^> i?*, when only the vertical pairs 
(I 3/2 — 2/1 I 3> I x 2 — x\ I ~ i?*) contribute. The scale i?* 
itself plays negligible role in the renormalization, since 
near the critical point in ID the total contribution comes 
from exponentially large range of distances. We thus can 
simply sew together the known RG descriptions Q (see 
also for the purely commensurate and the purely 

disordered cases: 



dK 
~d\ 



w 



dw 



{2-l/K)w (A < A*) 



(10) 



^§=Kw\ ^ = (3/2-l/K)w (A > A*) . (11) 

Here A = lni?, A» = lni?*; K = [tt-^/AJk] 1 (A s is the 
superfluid stiffness, k is the compressibility). The quan- 
tity [w(A)] 2 is proportional to the number of vortex pairs 
of the size ~ R in the area ~ R 2 of the xy-plane [^4) , 
and is much less than unity in the region of applicability 
of Eqs.(|l0|)-(|ll|). Strictly speaking, u>(A) experiences a 
jump at A = A*, but this is not important for our con- 
sideration. 

The critical value of A"(A = 00) for the destruction of 
superfluidity (divergency of AT(A) as A — > 00) is 1/2 in the 
commensurate case |p^| , |l3| ,p| and 2/3 in the disordered 
one Hence, if AT(A*) < 2/3, then it may turn out 
that K (00) < 2/3 also, and the system is in the superfluid 
phase though in the pure case it would be Mott insulator. 
In the weak-disorder limit one may neglect the difference 
between AT(oo) and AT(A*), provided the system is in SF 
phase [that is neglect the renormalization of K due to 
disorder]. The critical condition for the SF-BG transition 
then reads 



AT(A») = 2/3 . 



(12) 



To reveal the behavior of AT (A*) as a function of the 
Hamiltonian parameters consider the solution of the sys- 
tem (|Tc|) - near the critical value K = 1/2 from the 
Mott-insulator side: 

K = i + atan[4aA-7r/2] (\ K - - \ < 1) . (13) 

[A free constant in the argument of the tangent is fixed to 
be — 7r/2 (up to a small uncertainty on the order a <C 1) 
by the condition that at the microscopic scale A ~ 1 
the parameter AT(A) should be noticeably smaller than 
its mesoscopic value 1/2.] Near the critical point the 
quantity a depends on the parameters of the Hamiltonian 
of the pure system in the square-root fashion: 



oc vWpo - 1 



(14) 



Here p is any parameter of the Hamiltonian, the value 
p = Po corresponding to the point of SF-MI transition in 
the pure system (p > pq in the Mott phase). For exam- 
ple, one may choose p to be the strength of interparticle 
interaction. 

From Eq.Q it is seen that A"(A) does not differ essen- 
tially from 1/2 until A comes closely to the point n/Aa. 
At this point K diverges and we thus conclude that to a 
good accuracy 



A* — 



4a 



(15) 



A comment is in order here concerning the usage of 
Eq.(|l3|) at K = 2/3, i.e. out of the region of its formal 
applicability. Actually, at this value of K Eq.dTl) yields 
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only a rough order-of-magnitude estimate. This rough- 
ness, however, is essentially compensated by the sharp 
behavior of K(X) in the vicinity of the point A = ir/4a. 
So that the relative uncertainty of the value of A* is easily 
estimated to be on the order a <C 1. 

Combining Eqs. (Q), (|l4|), and ( |l5| ) we obtain the equa- 
tion for the phase boundary between superffuid and Bosc 
glass [W = Wsf-bg(p)]: 



Wsf-bg(p) ~ exp 



8by/p/pa - 1 



(16) 



where b is some constant. Eq.([l6]) demonstrates the ef- 
fect of disorder-stimulated superfluidity in the vicinity 
of the SF-MI transition of the pure commensurate sys- 
tem, which was observed previously in numerical studies 
lIlUPI . Note a very sharp behavior of the critical value of 
p for the SF-BG transition as a function of the strength 
of disorder. 

To provide the 'enhanced' proof of the existence of the 



BG in the weak disorder limit we compare Eq.(16) with 
the equation for the phase boundary between MI and 
a gapless phase (which subsequently will be identified 
with BG; until this is done the notations MI-BG and SF- 
BG should be understood as a convention). To obtain 
the latter equation one can take advantage of the exact 
relation taking place at the MI phase boundary (HI : 



A = W , 



(17) 



where A is the insulating gap of the pure system. This 
relation follows from the fact that the infinitesimal de- 
struction of the MI is due to the most favorable exchange 
of particles between exponentially rare Lifshitz's regions, 
where local configurations of disorder mimic uniform po- 
tential of the value W/2 or -W/2 An estimate for 
A one can obtain from the relation A ~ 1 /i?* || . So we 
have 



W-mi-bg (p) ~ exp - 



- 1 



and see that 



Wsf-bg (p) 
Wmi-bg (p) 



exp 



> 1 



(18) 



(19) 



Hence, in a weakly disordered commensurate ID system 
the strength of disorder sufficient to restore superfluidity 
is exponentially larger than that enough to turn Mott 
insulator into a gapless phase. This immediately proves 
the existence of BG between SF and MI, and identifies 
( |l6| ) and ( [L8| ) with corresponding phase boundaries. Ul- 
timately, this effect is due to the fact that A scales like 
1/i?*, while the critical amplitude of disorder supporting 
superfluidity is ~ 1 / V-R* ■ 



Let us now discuss recent results of numerical density- 
matrix RG study of ID commensurate disordered bosonic 
Hubbard model l|. In the weakly disordered case the 
authors observe a direct transition from SF to MI, in 
contrast to the case of a substantial disorder where BG 
intervenes between SF and MI, and to our asymptotically 
exact (at W — > 0) analytical results. 

It is known, however, that in the limit of weak disorder 
it is difficult to distinguish BG from MI by a finite-size 
scaling, because the difference is associated with expo- 
nentially rare Lifshitz's regions ||f[0]]- If the size of the 
system is not large enough, these regions are simply ab- 
sent with probability close to unity, and the scaling anal- 
ysis reveals nothing. [Maximal size of the system in Ref. 
H is of the order of 100 sites.] 

Therefore, one may suspect that the transition ob- 
served in Ref. at a small disorder is actually the SF- 
to-BG one. In this connection it is interesting to com- 
pare the numerical data with Eq. ( |l6| ) . Unfortunately, the 
phase diagram in Ref. || contains only few points in the 
weak-disorder region, too little to make possible an un- 
ambiguous fitting. However, one may check a non-trivial 
qualitative prediction of Eq.fll6|), the requirement that 
the strength of disorder on the SF-BG phase boundary 
be much larger than the Mott gap in the absence of dis- 
order. The data clearly demonstrate this feature. At any 
of the three presented points on the superffuid - insulator 
phase boundary in the weak-disorder region the strength 
of disorder is at least an order of magnitude larger than 
the insulating gap. 

We can allow also for the effect of a small deviation 
from commensurability. To this end we should add to 
the right-hand side of Eq.(|4|) the term 2ir8'n/ \xi — X2), 
where Sn^ is the deviation of the density from the com- 
mensurate value. This term is negligible for a vortex pair 
of a radius much less than 



Rr 



(20) 



(8v is the deviation of the filling factor), and leads to 
the cancellation of contributions of pairs with radii R c < 
R < i?*, provided R c < i?*. At R c > R* this term 
plays no role at all. Hence, at R c < R*, we just need to 
replace A* by lni? e in Eq. (|io|) and set if (A*) = K{\aR c ) 
as the initial condition for Eq. ([ll]) . Consequently, for the 
SF-BG phase boundary we obtain the equation 



exp 



\W\\8v\) . (21) 



Note that in the variables p and 8v (at fixed W) the 
phase boundary has an independent of W cusp-like shape 
except for the tip of the cusp cut off by disorder. 

A remark is in order here concerning fractional filling 
factors. The fact that the point of SF-MI transition of 
a pure commensurate system turns out to be the limit- 
ing point for both SF-BG and BG-MI phase boundaries 
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of the disordered system is characteristic of only integer 
filling factors, because only in this case the critical value 
of the parameter K for the pure system, K = 1/2, is less 
than that of disordered system, K — 2/3. In the vicin- 
ity of SF-MI transition at a fractional filling, where the 
transition occurs at K — q 2 /2 (q is the denominator of 
the filling factor) |l5|,|l3|], an infinitesimal disorder takes 
the system away from the superfluid region. 

Summarizing, we have presented a description of the 
superfluid - Bose-glass phase boundary in ID weakly 
disordered system near the point of superfluid - Mott- 
insulator transition of the pure system with integer fill- 
ing factor. The results demonstrate the effect of stim- 
ulation of superfluidity by disorder, previously observed 
numerically. In the limit of weak disorder the strength 
of disorder sufficient to restore superfluidity from Mott 
insulator turns out to be much larger than that enough 
to turn Mott insulator into Bose glass. This provides an 
explicit proof of the fact that at arbitrarily small disorder 
the Bose glass always intervene between superfluid and 
Mott phases. 
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